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Self similarity of two point correlations 
in wall bounded turbulent flows 

By J. C. R. HUNT 1 , P. MOIN 2 3 , R. D. MOSER 3 AND P. R. SPALART 3 


Computational results of Moin & Moser (1987) and Spalart (1986) for two-point 
correlations of the normal v component of turbulence at two points y, 1/1 (yj > 
y), from the rigid walls bounding turbulent channel and boundary layer flows for 
Reynolds numbers 3200 and 7000 are shown to have an approximately self-similar 
form, when plotted in terms of yjy\. It is found that 

HvHyi) = ,/j/v 

v 2 (yi) V 

where 0 < y/yi <l,/(0) = 0, /(1) = 1, and where / is approximately independent 
of 3 /i , for r/i ranging from about 20 i//u x to half the channel width; 

/ % 2(y/yi) 2 - (y/j/i) 3 ± .1. 

The same kind of self similarity has been predicted for and measured in shear 
free boundary layers. But in that case, where / % y/y l5 the mechanism is one of 
‘blocking’ or ‘splatting’ at the wall. In these sheared wall layers, the shear also 
has an important effect. There are important implications from this research for 
modeling wall bounded shear flow. 


1. Introduction and Objective 

The structure of turbulence at a height y from a wall is affected by the local mean 
shear at y, (§^), by the direct effect of the wall on the eddies, and by the action 
of other eddies close to or far from the wall. Some researchers believe that a single 
one of these mechanisms is dominant, while others believe that, these effects have 
to be considered together. 

It is important to understand the relative importance of these effects in order to 
develop closure models, for example for the dissipation or for the Reynolds stress 
equation, and to understand the eddy structure of cross correlation functions and 
other measures. The specific objective of this research project was to examine the 
two point correlation R vv of the normal velocity component v near the wall in 
a turbulent channel flow and in a turbulent boundary layer. This component of 

1 Univ of Cambridge 

2 Stanford University 

3 NASA Ames Research Center 


PRECEDING PAGE BLANK NOT FILMED 



26 


J. C. R. Hunt, P. Moin, R. D. Moser and P. R. Spalart 


V L + V L “ V/Vl V L 
(COMBINED EFFECT 



(a) (b) 


FIGURE 1. Diagram to illustrate eddy motion near the ground, (a) typical struc- 
ture of an updraft or “thermal” plus the “image” updraft, ( b ) the relation between 
the velocity at t/j and y. The large eddy L , with velocity vl, is centered at height 
yj. The small eddy S , with velocity (’5 is centered at height y. (Their “images” 
are V and S') The profiles are shown of the vertical velocity of the large and small 
eddies, of their images and of the combined effect of both. 


turbulence is the most sensitive to the relative effects of shear (which amplifies v) 
and the blocking effect of the surface (where v = 0 , even in inviscid flow). 

Recent research on shear free turbulent boundary layers, (such as occur in thermal 
convection between boundaries or in turbulence near a free surface or turbulence 
near a density inversion layer) has shown how the blocking effect leads to a self- 
similar form for R vv when expressed as a ratio with u 2 (yi) (i.e., normalized at the 
upper point), 


= A* ) - i »<» 

r 2 (yi ) Vi 2/i 


The theory for the SFBL is valid when yi is much less than the turbulence scale 
far from the boundary. The explanation is given with the aid of figure 1 . Let there 
be a large eddy (L) centered with maximum velocity at y\ and a small eddy centered 
at y. Then the velocity at yj, v(yj ) % v L . The effect of the small eddy at y is small 
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FIGURE 2. Measurements in the convective atmospheric boundary layer of the 
cross correlation of v at heights y and t/j , normalized by v 2 at j/j . The results are 
compared to the approximate form of the theoretical predictions of Hunt (1984), 
Rw ^ y/yi' (~Zi/L indicates the state of the convective boundary layer, the higher 
— Zi/L, the stronger the convection and the weaker the shear.) From Hunt, Kainal 
& Gaynor, (1987). 


if y is less than about t/j/2. However, the velocity at y 

v(y) = v, + (y/yi )v L (1.2) 

has two components, from the small eddy and also from the large eddy. Because 
the vertical dependence of the large eddy is blocked by the surface, this component 
is reduced by a factor of about (y/yi) for high Reynolds number turbulence. One 
can imagine an image vortex underneath the surface. Since the correlation between 
v s and Vi is small if y < t/i/2 the correlation between v(y) and v(yi) and thence 
the correlation, normalized at the upper point (N.U.P.) is 

Rw ~ y/yi. ( 1 - 3 ) 

It is interesting that the theory seems to agree with measurements for the atmo- 
spheric boundary layers during thermal convection. See figure 2 from Hunt, Kaimal 
& Gaynor 1987. The same general idea might be appropriate for a wall bounded 
shear flow at moderate Reynolds number, but now the velocity at height y is not a 
simple function of y/y^. In general we expect the component of vi to be given by 
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a function g[{y - yi)/L v (yi)], depending on the distance between y and t/] and the 
scale of turbulence L v (yi) at yj. Therefore, we might predict: 

Rv»*g((y-yi)/L v ( yi )) (1.4) 

But in a wall-bounded flows L (,,) (yi) is proportional to yj, (at least for unstable 
and neutral, but not for stably stratified flows). Therefore, we might expect, 


9{ 


y - m 

L v (y ,) 


) =d{ 


y - y i 

ay i 


) 


where L v (y \ ) « ayi, which can be written as another function, i.e. 

Rw ~ /(y/yi). 


(1.5) 


This is the first hypothesis to be tested. 

We have argued that the wall “blocks” the normal component of the large-scale 
eddies centered above the wall at y = yj. In a shear flow the streamwise u and 
normal v components of the turbulence are correlated at the height yj, i. e., R uv = 
uv(y 1 ) ^ 0. Therefore the cross-correlation between the normal velocity v(y) at y 
(< y \ ) and the horizontal velocity at yi, tt(yi), should steadily decrease near the 
wall as y/yj — ► 0. If the scaling argument of (1.5) is valid one should expect that, 
in the log layer, where y << yj, the u — v correlation normalized by uv at y t has 
the form: 


R 


UV 


Ruviyi i y ) ^ r ( y_\ 
Ruv(yi,yi) 


( 1 - 6 ) 


This is the second hypothesis to be tested. 

In general, correlations involving the horizontal component are affected by the 
inactive or irrotational motions. Consequently, the presence of the wall exerts a 
weaker influence on these correlations. 

The computations of the structure of homogeneous t urbulence in a uniform, mean 
velocity gradient have shown that the main effect of the mean shear is to reduce the 
scale of the turbulence in the spanwise or 2 direction (Townsend 1976). Unlike many 
statistical effects this one is so strong that this channeling of turbulent eddies can be 
seen in the instantaneous pattern of streamlines, as derived from flow-visualization 
studies and direct simulations (e. g., Lee, Kim and Moin 1987). 

The picture of these eddies in shear flow, indicated in figure 3, looks very different 
from the conventional circular vortex-like eddy of homogeneous isotropic turbulence; 
such as illustrated in figure 1. One suggestion is to represent the eddies as vertically- 
elongated structures with a defined spanwise scale a 3 . 

This suggests that the structure of the large eddy is approximately described by 


Self Similarity of Two Point Correlations 


29 




FIGURE 3. Schematic of spanwise structure of eddies in the boundary layer, 
showing how the spanwise scale (a 3 ) of the large eddy ( L ) at y 3 determines the 
spanwise scale closer to the wall (at y ). Also, if y 3 is smaller, a 3 is also smaller. 


(as opposed to the isotropic eddy structure where v % /[(rj + r\ )/L 2 ]vi). Here 
r 2 = (j/i — y) and r 3 are the distances from the center of the eddy. Thus at a height 
y, and spanwise displacement z , 

v(y,r 3 ) = v, + v L f v (—) g( — ) 

2/i ® 3 

and therefore 

v(y,r 3 )v(y u Q) % 

This model also implies that 


Rw{y^ ) 


Rvv{y,r 3 ) _ ^(y,r 3 )v(yi,0) _ ^ , 

Rw(y,Q) v(y,0)v(j/i,0) 9 «3 


( 1 . 8 ) 


where the function g is independent of 2 . If the hypothesis (1.6) is valid we would 
also expect that 


Rvu{y,r 3 ) 


R^ui^n) 

Rvu(y,0) 


t’(y,r 3 M 2/i,0) ^ rs 
u(j/,0)u(yi,0) 9 «3 


Since the eddy structure (1.7) occurs only on strong shear flows, it is natural to 
suppose that a 3 depends on the mean velocity gradient dU/dy and the vertical 
turbulence intensity v' in the log layer. Near the wall the lateral structure is likely 
to be determined by instabilities within the wall layer. So we postulate that 


<*3(3/1) ~ <* 3 


dU / dy(yi ) 


+ ct w 



(1.9) 


This is the third hypothesis to be tested here. 
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FIGURE 4. Cross correlation of v at heights y and yi , normalized by v 2 at yi 
computed from direct numerical simulations of the zero pressure gradient boundary 
layer (Spalart 1987), and the plane channel (Moin k Moser 1987). Also shown are 
the theoretical predictions of Hunt (1984) ( R vv ~ y/yi )• 


2. Preliminary results 

In figure 4 we present a graph of R V v(y/y l) combining the results of the compu- 
tations of the zero pressure gradient (Z. P. G.) boundary layers (Spalart 1986) and 
of the channel flow (C.F.) (Moin k Moser 1987). The range of data is as follows: 
ZPG: Sur/u = 300 (where 8 is the boundary layer thickness) and 

y + = 10, 22, 40, 97, 190 

Channel Flow: 8u r /v — 180 (where 8 is the channel half width) and 


y-\u T /v = 33, 59, ... 180 


Self Similarity of Two Point Correlations 


31 


In figure 4 the results for the channel alone are plotted including the exceptionally 
small value of t/j —hvju r . 

It appears that for these two wall bounded shear flows, the self similar plot of 
R vv is a good description of the measurements. Comparing figures 2 and 4 indicates 
as good a ’collapse’ as observed for shear free boundary layer. (But remember that 
figure 2 is a plot of experimental points in the atmosphere!). Note that the similarity 
hypothesis is more accurate for smaller values of y, as expected since the assumption 
that the small eddies at y and the large eddies as y j are uncorrelated, (i.e. vJvl ~ 0) 
is more valid when yjy\ is small. 

It is particularly surprising that approximately the same curve describes the dis- 
tribution of R vv for points both within and well above the viscous sublayer. However 
very close to the surface within the sublayer we must expect that, since v oc y 2 as 
yu T / v -*• 0, R vv — > (y/yi) 2 as yu r /v -» 0. 

However the computation of Kim, Moin or Moser (1987) have shown that some 
vertical eddying motion exists on a scale even smaller than hv /u r (because v is not 
exactly proportional to y 2 ). This is quite consistent with the fact that the two point 
correlation R vv is greater than (y/yi) 2 when y^u T jv = 5. 

These results show that there is a significant difference in the measured value 
of R vv between these shear boundary layers (figures 4 and 5a) and the shear-free 
boundary layers in figure 2. They show that the effect of shear is to reduce the cor- 
relation length of the normal velocity in the normal direction. (But it is important 
to note that the smallest scale of v in a shear flow is in the spanwise or z-direction, 
Townsend (1976). So these curves of R vv do not give a basis for estimating dissipa- 
tion or the dissipation length scale.) 

Figure 5 b is added to show that if the conventional two point correlation is plotted 
against y/yj, the points do not tend to zero as y/yj — * 0 and the curves do not 
have any general pattern. 

In figure 5c we present the cross correlation of the Reynolds stress R uv as defined 
by (1.6). These curves for different values of ( y\/8 ) show that R uv is not a universal 
function of (y/yi) over the whole channel width. However, the three sets of curves 
for which yj is in the log layer, i.e., 30 < y+ < 100 do exhibit a strong degree of 
similarity- (recall that the abscissa is being stretched by a factor of 3, and yet the 
curves are self similar). 

It is interesting to note that the shape of the self similar curves of R uv is markedly 
different to the curves for the normal velocity correlation, R vv . The correlation is 
higher. 

In figure 6, we present the cross correlation for the vertical velocity separated by 
a normal and spanwise spacing normalized by the correlation at the same height, 
but zero spanwise spacing i.e., 


R (y^3,yi) _ t’(y,r 3 )u(y 1 ,0) 
R vv (y,0,yi) v(y,0)v(yi,0) 


We also plotted the correlation of u and v defined as 


Rvx,(y,r 3 ;yi) = R„ u (y,r- 3 ;y 1 )/^ t , v (y,0,y 1 ) 
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FIGURE 5. Normalized cross correlations of velocities at y and t/i from tur- 
bulent plane channel flow (Moin & Moser, 1987) for various values of y\. (a) 

Hy)vM/v 2 (yi), (b) v(y),v(y 1 )/yjv 2 (y)v 2 (y 1 ), (c) u{y)v{y 1 )/u(y 1 )v(y i ). 


Self Similarity of Two Point Correlations 


33 


© 







- 1.0 


JJb) 

0 


r 3 /a 3 


— r- 

2 


3 



FIGURE 6. Normalized cross correlations of velocities at y and yj with separation 
in 2 (r 3 ) fr om turbulent pl a ne channel flow (M oin fe Moser, 1987) for various v alues 
oft/,, (a) v(y,r 3 ),v{y u 0)/v(y,0)v(y u 0), (b) v{y,r 3 ), u( yi ,0)/v(y,0)u( yi ,0). The 
value of a 3 is chosen for each curve to be the r 3 location at which the curve passes 
through 1/e. 
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FIGURE 7. Values of «3 as a function of y/yt and yj, (a) as used to scale R vv 
(figure 6a), (b)) as used to scale R vv (figure 66). 


The curves are plotted as functions of r 3 /a 3 for different, y and y \ , where <13 is 

defined as the value of r 3 where R vv = 1/e = 0.36 The results show firstly, the 
small variation of the form of the spanwise structure and the negligible variation 
of the spanwise scale of the eddies in these boundary layers. Using this particular 
correlation emphasizes this point quite nicely. This approximate invariance is found 
for values oiyi/6 < 0.8. It is not true for yj at the centerline, yi/6 23 1.0. 

Secondly, these results show how the scale a 3 increases with yj. In Figure 7, we 
have plotted against y"*" . It appears that this scale is of the order of 9 wall units 
near the wall and then begins to increase when y + > 10. This would be consistent 
with the ideas suggested in the introduction. A satisfactory curve fit could be (figure 
8 ) 


dU + 

( L4 dy+ +7 ^ 

0.3y + in the log layer. 


flS 
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FIGURE 8. Values of a 3 as used to scale R vv (figures 6a and 7a), as a function of 
the reference height y 3 . x - simulation results, o - the model a 3 = 7 + 


This value of is of the same order but a little greater than the dissipation length 
scale L ( (based on v 2 in the log layer) where 

L ( = ~ 0.18y 

v 2 

Other results have been computed for R UV ,R UU . They also show equally strong 
channeling of the spanwise structure; though the value of a 3 for R uu is about twice 
as great as for R vv . 

3. Implication and further work 

The preliminary results show that even in the inhomogeneous turbulent boundary 
layer, the two-point correlation function may have self similar forms. The nature 
of these self similar functions can be inferred by using rapid-distortion theory. The 
results shows that the effects of shear and of blocking are equally important in 
the form of correlation functions for spacing normal to the wall. But for spanwise 
spacing, we have found that the eddy structure is quite different in these shear 
flows; this aspect of eddy structure is largely controlled by the shear and perhaps 
by small scale structures very close to the wall. So any theory for the turbulent 
structure must take both these effects into account. 

The results suggest further study 
. Comparison with laboratory and atmospheric measurement. 

. The effects of curvature and pressure gradient should be investigated. We would 
still expect to see these self similar two point correlations, 
c. Further theoretical calculations should be done using RDT Theory for uniform 
shear near a wall. 
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d. The two point correlation functions are likely to be self similar for other compo- 
nents and in other direction (eg. for the spanwise, z, spacing). 
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